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CN ■ Abstract 
(N 

We present a search for CO (3-2) emission in SDF-26821, a BzK -selected star-forming galaxy (sBzK) at z = 
2.044, using the 45-m telescope of the Nobeyama Radio Observatory and the Nobeyama Millimeter Array. We do 
not detect significant emission and derive 2 a limits: the CO luminosity of L' co < 3.1 x 10 10 K km s _1 pc~ 2 , the 
, ratio of far-infrared luminosity to CO luminosity of Lfir/-^co - > ^7 Lq (K km s" 1 pc~ 2 ) -1 , and the molecular 

gas mass of Mh 2 < 2.5 x 10 10 Mq, assuming a velocity width of 200 km s" 1 and a CO-to-H2 conversion factor 
of «co = 0.8 Mq (K km s" 1 pc~ 2 ) -1 . The ratio of Lfir/ L' co , a measure of star formation efficiency (SFE), is 
comparable to or higher than the two z ~ 1.5 sBzKs detected in CO (2-1) previously, suggesting that sBzKs can 
have a wide range of SFEs. Comparisons of far-infrared luminosity, gas mass, and stellar mass among the sBzKs 
(— > , suggest that SDF-26821 is at an earlier stage of forming stars with a similar SFE and/or more efficiently forming 

stars than the two z ~ 1.5 sBzKs. The higher SFEs and specific star formation rates of the sBzKs compared to local 
spirals are indicative of the difference in star formation modes between these systems, suggesting that sBzKs are not 
just scaled-up versions of local spirals. 
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, 1. Introduction et al. 2004). BzA-selected galaxies (BzKs) are classified into 

two categories, star-forming galaxies (sBzKs) and passively- 
Molecular gas, fuel for star formation, is a fundamental com- evolving galaxies (pBzKs), based on their colors. CO obser- 
■ ponent of galaxies. Observations of molecular gas, traced vations of BzKs and BXs have been performed by Daddi et 
| by CO emission lines, provide important information on the al. (2008a) (hereafter D08) and Tacconi et al. (2008). D08 
evolution of galaxies: e.g., star formation efficiency and the successfully detected CO (2-1) emission from two sBzKs at 
timescale and the formation phase of galaxies. In the forma- z ~ 1.5, yielding molecular gas mass of Mh 2 ~2x 10 10 Mq. 
tion history of galaxies in the universe, the era of z ~ 1-3 is The large amount of molecular gas with relatively low stellar 
thought to be a turning point, since the cosmic star formation mass (M* = 5 x 10 10 and 8 x 10 10 Mq, adopting a Chabrier 
rates peaked and the majority of stellar mass had formed by (2003) initial mass function (IMF)) and high star formation 
z ~ 1 (e.g., Dickinson et al. 2003). Therefore it is important to rates (SFRpm = 130 and 235 Mq yr -1 ) imply that gas-rich 
reveal the gas environment of this era to probe the star forma- galaxies in the process of formation prevail in the early uni- 
tion activities in galaxies. However, the detection of CO line verse, considering the fact that the volume density of sBzKs is 
at this era has been made mainly in apparently bright galaxies a factor of 10 larger than that of SMGs (e.g., Daddi et al. 2004; 
such as radio galaxies, QSOs, and bright submillimeter galax- Chapman et al. 2005). On the other hand, Tacconi et al. (2008) 
ies (SMGs), or in lensed galaxies (see Solomon & Vanden Bout did not detect CO emission in a BzK and two BXs at z ~ 2, pro- 
2005 for a review). For a comprehensive understanding of the viding upper limits of molecular gas mass of 5-7 x 10 9 Mq (2 
galaxy evolution, it is necessary to observe molecular gas in a). This is inconsistent with galaxies in the local universe be- 
galaxies which dominate the galaxy population in the universe, cause we expect a gas mass of a few 10 10 M© for a galaxy with 
The BzK color selection technique using two colors, B — z' SFR of ~ 200-300 Mq yr _1 assuming a standard Schmidt law 
and z' — K, provides a A"-limited (i.e., mass-limited) sample (Schmidt 1959; Kennicutt 1998). 

of galaxies at 1.4 < z < 2.5 (Daddi et al. 2004). This tech- This may be caused by differences in physical properties 

nique is less affected by dust extinction and effectively selects of these systems such as CO-to-H2 conversion factors, excita- 

dusty starburst galaxies compared to samples of z ~ 2 UV- tion conditions of CO, relation between gas surface density and 

selected star-forming galaxies, BX and BM galaxies (Steidel SFR, and star formation modes. It is essential to increase the 
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number of CO observations of these high-redshift star-forming 
galaxies to examine the physical properties. 

In this letter, we report a search for CO (3-2) line emission in 
an sBzK detected in the Subaru Deep Field (SDF; Kashikawa 
et al. 2004). Our target, SDF-26821, is selected from deep 
if-band imaging (Hayashi et al. 2007) and confirmed to be at 
z = 2.044 ± 0.002 by the Subaru/MOIRCS near-infrared spec- 
troscopy (Hayashi et al. 2009, hereafter H08). The FIR lumi- 
nosity is estimated to be Lpm, = 1-7 x 10 12 L© from the 24 fjm 
flux obtained with the Spitzer/MlPS observations (Ly et al. in 
prep.), using a template SED of local galaxies (Chary & Elbaz 
2001), with about factor of a few uncertainty (Papovich et al. 
2007). Thus, this galaxy is classified as an Ultra Luminous 
Infrared Galaxy (ULIRG). The SFR is the largest among the 
spectroscopic sample: 370, 730, and 300 M@ yr _1 derived 
from rest-frame 1500 A corrected for dust extinction using the 
Calzetti extinction law, extinction-corrected Ha, and FIR lu- 
minosities, respectively (H08). Although there is a possible 
contribution from an active galactic nucleus (AGN), the con- 
tribution appears to be insignificant (H08). The stellar mass of 
3.4 x 10 10 M o is derived from the SED fitting (with about 30% 
of fitting error) assuming a Salpeter (1955) IMF with lower and 
upper mass cutoffs of 0.1 M Q and 100 Mq and a continuous 
burst. The large SFR compared to the relatively small stellar 
mass suggests that SDF-26821 is in its forming phase and may 
contain a large amount of molecular gas. 

In §2 we describe our CO (3-2) observations and data re- 
duction. In §3 we present the results and physical quantities 
obtained. In §4 we compare CO luminosity, FIR luminosity, 
and stellar mass of SDF-26821, sBzKs, and other galaxy pop- 
ulations and discuss the difference in star formation. Where 
necessary, we adopt a cosmological parameter of Hq = 70 km 
s- 1 Mpc" 1 , Q M = 0.3, and Q A = 0.7. 

2. Observations and Data Reduction 

CO (3-2) observations were initially conducted with the 45- 
m telescope of the Nobeyama Radio Observatory (NRO) on 
March 13-15, 2008. The observation frequency was set to the 
redshifted CO (3-2) line of 113.599 GHz at the upper side- 
band of the T100 receiver, a waveguide-type dual -polarization 
sideband-separating SIS receiver (Nakajima et al. 2008). The 
S100 receiver was also used for gain calibrations of the T100 
receiver. The half-power beam width (HPB W) of the telescope 
for two polarizations were 15". 5 ± 0".l and 17". 6 ± 0".3, and 
the main-beam efficiency was 34% and 35%. We used the 
AOS-Wide spectrometer with 2048 channels and a bandwidth 
of 250 MHz. The pointing was checked every hour using SiO 
masers. The sky emission was subtracted by position switch- 
ing. The intensity calibration was performed by the chopper 
wheel method, and the absolute intensity calibration was per- 
formed by observing IRC+10216. The system temperature 
was in the range of 300-500 K in SSB. 

We then conducted observations with the Nobeyama 
Millimeter Array (NMA) on May 19-27, 2008. We used the 
most compact array configuration, D-configuration (baseline 
length: 15-82 m), with five 10-m antennas. The Ultra Wide 
Band Correlator (UWBC; Okumura et al. 2000) was used with 
128 channels and a bandwidth of 1024 MHz, covering ~ 2700 




Right Ascension (J2000) 

Fig. 1. NMA contour map around the position of SDF-26821 (cross) 
at 113.63 GHz with a velocity resolution of 84.3 km s" 1 (1 channel). 
Contours are —2 a, —1 a, 1 a, and 2cr (1 a = 5.2 mJy). The large circle 
represents the field of view (62" diameter). The synthesized beam is 
shown on the bottom left. 
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Fig. 2. NMA spectrum at the position of SDF-26821 with a velocity 
resolution of 84.3 km s~ 1 . The vertical solid and dashed lines indi- 
cate the frequency range of redshifted CO (3—2) line expected from the 
near-infrared spectroscopy. 

km s -1 . The bandwidth was wide enough even allowing for 
the uncertainty of observation frequency (~ 150 MHz) due to 
the uncertainty of redshift. The observation frequency was set 
to the redshifted CO (3-2) line of 113.599 GHz at the upper 
sideband, 200 MHz shifted downward from the center of the 
band to prevent spurious features. The field of view is 62" 
in diameter. The QSO 1308+326 was observed every 20 min- 
utes for visibility calibrations, and 3C 273 were observed for 
bandpass calibrations. Data reduction was carried out using 
the UVPROC2 software developed at the NRO (Tsutsumi et al. 
1997) and the AIPS package developed at the National Radio 
Astronomy Observatory (NRAO). The total on-source integra- 
tion time was 16.4 hours. The synthesized beam size was 
5". 4 x 4". 6 (position angle of —44 ) with natural weighting. 

3. Results 

No significant CO line and continuum emissions are de- 
tected with either observations. The contour map and spec- 
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trum obtained with the NMA are shown in Figures 1 and 2. 
The rms noise levels for the 45-m and NMA observations are 
1 1 mJy and 6.8 mJy (50 km s _1 resolution), respectively. The 
rms noise level of continuum emission with NMA is 1 . 1 mJy. 
This is more than an order of magnitude higher than expected 
from the SED model and SFRs, and does not impose tight con- 
straints on the SED. The noise level for NMA is measured at 
the source free region of the contour map avoiding the source 
signal expected at the phase reference center. In what follows, 
we discuss 2 a limits using the rms noise level obtained with 
the NMA observations where deeper sensitivity was achieved. 
Following Solomon et al. (1992), the CO luminosity is 
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where L' co is measured in K km s 



Sqo is the observed 
CO flux in Jy, Av is a velocity width in km s -1 , and Dl is 
the luminosity distance in Mpc. Assuming a velocity width of 
200 km s" 1 based on the CO spectra of BzKs reported by D08, 
the 2 a upper limit of CO (3-2) luminosity is L' co < 3.1 x 10 10 
K km s _1 pc~ 2 . 

The molecular gas mass is given by Mh 2 = acoL' co , wnere 
aco is a CO-to-H2 conversion factor. We assume that the gas 
is optically thick and thermalized, and a CO (3-2)/CO (1-0) 
luminosity ratio of unity. By adopting a conversion factor of 
a co = 0.8 M (K km pc 2 )" 1 , the standard value for ULIRGs 
(Downes & Solomon 1998), the 2 a upper limit of molecular 
gas mass is Mh 2 < 2.5 x 10 10 M . This is lower than the 
median value of (3.0 ± 1.6) x 1O 1O M obtained in the sample 
of 12 SMGs (Greve et al. 2005). 

4. Discussion 

The CO luminosity and the FIR luminosity are measures of 
molecular gas mass and SFR, respectively, and therefore the 
ratio of Lpm/ico indicates how efficiently stars are formed 
from molecular gas. It is also used as an indicator of star 
formation efficiency (SFE; Young et al. 1986). Figure 3 
shows Lfir/^co ra ti°s as a function of redshift using vari- 
ous populations of galaxies; local spiral galaxies, local lumi- 
nous infrared galaxies (LIRGs, L FIR = lO n ~ 12 L ), ULIRGs 
(£fir > 10 12 £©), sBzKs, Lyman break galaxies (LBGs), 
SMGs, QSOs, and high-redshift radio galaxies (HzRGs). The 
2 a lower limit of SDF-26821 is L FIR /L' co > 57 L (K km 
s _1 pc 2 ) -1 . Considering it is a lower limit, SDF-26821 has 
a SFE comparable to or higher than those of two sBzKs of 
D08 (U±H and50±^ L (Kkms" 1 pc 2 )" 1 ), suggesting that 
sBzKs could span a wide range of SFEs. The ratios of SDF- 
26821 and D08 sBzKs are higher than the average value for 
local spirals (26 ± 3 L (K km s" 1 pc 2 ) -1 ), suggesting that 
sBzKs are not just scaled-up versions of local spirals. 

In Figure 4, we plot the ratio of FIR luminosity to molecular 
gas mass (Lftr/Mh 2 , another form of SFE) versus the ratio 
of stellar mass to molecular gas mass (M*/Mh 2 ) of ULIRGs, 
sBzKs, LBGs, and SMGs. For comparison, we show M 51 
as a typical spiral, M 82 as a typical starburst, and the Milky 
Way. The 2 a lower limit of SDF-26821 is represented by 
an arrow. In this SFE-M*/Mh 2 diagram, the ratio of verti- 
cal axis to horizontal axis provides specific star formation rates 
(SSFRs; SFRs per unit stellar mass) by using a conversion from 
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Fig. 3. Star formation efficiency (Lyir/Lqq) as a function of red- 
shift. The 2 a lower limit of SDF-26821 is represented by an up- 
ward arrow. For comparison, ratios of various galaxy populations are 
plotted: local spirals (Young et al. 1996), LIRGs, ULIRGs (Sanders 
et al. 1991; Solomon et al. 1997), BzKs (D08), LBGs (Baker et al. 
2004; Coppin et al. 2007), SMGs (Greve et al. 2005; Daddi et al. 
2008b; Tacconi et al. 2008), QSOs, and HzRGs (Solomon & Vanden 
Bout 2005; Carilli et al. 2007; Maiolino et al. 2007; Coppin et al. 2008, 
and references therein). 

FIR luminosity to SFR (Kennicutt 1998). The SSFR is thought 
to be an indicator of current star-forming activity, and the in- 
verse of the SSFR is related to the mass doubling time. The 
SSFR of SDF-26821 (8.8 Gyr -1 ) is a factor of six higher than 
those of D08 sBzKs (1.4 and 1.6 Gyr -1 ), suggesting that the 
mass doubling time of SDF-26821 is shorter than those of D08 
sBzKs. Compared to the other galaxy populations, SDF-26821 
and D08 sBzKs have similar SSFRs to the average for the sam- 
ple of SMGs (13 ± 4 Gyr" 1 ) and ULIRGs (2.4 ± 0.8 Gyr" 1 ), 
respectively, and more than an order of magnitude higher than 
those of M51 and the Milky Way. Taking into account the 
information of molecular gas mass, the location of a galaxy 
within this diagram shows the states of star-formation in galax- 
ies (SFE, SSFR, and the fraction of stellar mass to molecular 
gas mass). Figure 4 suggests that SDF-26821 is a different 
system from D08 sBzKs: (i) if the actual CO luminosity of 
SDF-26821 is close to the 2 a upper limit, then the SFE is 
comparable to those of D08 sBzKs. On the other hand, the 
fraction of stellar mass to molecular gas mass of SDF-26821 is 
a factor of five smaller than those of D08 sBzKs and an order 
of magnitude smaller than M5 1 and the Milky Way. The high 
fraction of molecular gas mass to stellar mass in SDF-26821 
could relate to an earlier phase of forming stars compared to 
D08 sBzKs (e.g., Frayer & Brown 1997). This is supported 
by the inferred young age of 50 Myr derived from the SED 
model fitting (Motohara et al. in prep.), (ii) If the actual CO 
luminosity of SDF-26821 is smaller than the 2 a upper limit, 
then the SFE is higher than those of D08 sBzKs and closer to 
SMGs. This suggests that SDF-26821 is similar to massively 
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Fig. 4. The ratio of FIR luminosity to molecular gas mass versus the 
ratio of stellar mass to molecular using samples of distant galaxies 
(sBzKs, LBGs, and SMGs) and local galaxies (the Milky Way, M 51, 
M 82, and ULIRGs). The 2 cr lower limit of SDF-26821 is represented 
by an arrow. The Lfir/-Mh 2 ratio is a measure of star formation ef- 
ficiency (SFE). Diagonals show constant specific star formation rates 
(SSFRs; SFRs per unit stellar mass). SFRs are derived from FIR lumi- 
nosities (Kennicutt 1998) and stellar masses are converted to a Salpeter 
IMF. Note that we use dynamical mass for ULIRGs and baryonic mass 
for M 5 1 and M 82 instead of stellar mass, and they give upper limits on 
stellar mass. We use CO-to-H2 conversion factors of aco = 4.6 Mq 
(Kkmpc 2 ) -1 (Galactic value; Solomon & Barrett 1991)forM51 and 
M 82, and o c o = 0.8 M G (K km pc 2 )" 1 for ULIRGs and distant 
galaxies. The ULIRG samples are taken from Tacconi et al. (2002), 
SMG samples are taken from Tacconi et al. (2008) and Daddi et al. 
(2008b), and sBzKs and LBGs are the same as in Figure 3. The data of 
FIR luminosities and molecular gas masses are obtained from Chapman 
et al. (2005) for SMGs of Tacconi et al. (2008), Solomon et al. (1997) 
for ULIRGs, Young et al. (1996) for M 51 and M 82, and Scoville & 
Good (1989) and Dame (1993) for the Milky Way. The stellar masses 
are taken from Silva et al. (1998) for M 51 and M 82, and Flynn et al. 
(2006) for the Milky Way. 

star-forming galaxies at high redshifts and could be a 'scaled- 
down' version (i.e., lower gas and stellar mass) of SMGs . 

Our CO observation reveals that sBzKs could have a wide 
range of SFEs. It is possible that sBzKs are not just scaled- 
up versions of local spirals and are more like vigorously star- 
forming galaxies at high redshifts. In order to examine the 
evolution process of these high-redshift star-forming galaxies, 
more observations of molecular gas are needed. However, the 
sensitivity of current instruments is restricted to only bright 
sources. ALMA will enable us to study the gas environment 
of normal galaxies at high redshift and to examine the galaxy 
evolution in terms of molecular gas, which has been less un- 
derstood thus far. 

We would like to acknowledge Takeshi Sakai, Taku 
Nakajima, and the members of the NRO for observational sup- 
port. We are also grateful to Emanuele Daddi for kindly provid- 
ing detailed data for the BzKs. This work has been supported 
in part by the Grant-in-Aid for Scientific Research (19540238) 



from Japan Society for the Promotion of Science (JSPS). B. H. 
and M. H. are financially supported by a Research Fellowship 
from the JSPS for Young Scientists. 

References 

Baker, A. J., Tacconi, L. J., Genzel, R., Lehnert, M. D„ & Lutz, D. 

2004, ApJ, 604, 125 
Carilli, C. L., et al. 2007, ApJL, 666, L9 
Chabrier, G. 2003, ApJL, 586, L133 

Chapman, S. C, Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 
622, 772 

Chary, R, & Elbaz, D. 2001, ApJ, 556, 562 

Coppin, K. E. K., et al. 2007, ApJ, 665, 936 

Coppin, K. E. K., et al. 2008, MNRAS, 389, 45 

Daddi, E., Dannerbauer, H., Elbaz, D., Dickinson, M., Morrison, G., 

Stern, D., & Ravindranath, S. 2008, ApJL, 673, L21 (D08) 
Daddi, E., Cimatti, A., Renzini, A., Fontana, A., Mignoli, M., 

Pozzetti, L., Tozzi, P., & Zamorani, G. 2004, ApJ, 617, 746 
Daddi, E., et al. 2008, arXiv:08 10.3 108 
Dame, T. M. 1993, Back to the Galaxy, 278, 267 
Dickinson, M., Papovich, C, Ferguson, H. C, & Budavari, T. 2003, 

ApJ, 587, 25 

Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 

Flynn, C, Holmberg, J., Portinari, L., Fuchs, B., & JahreiB, H. 2006, 

MNRAS, 372, 1 149 
Frayer, D. T, & Brown, R. L. 1997, ApJS, 113, 221 
Greve, T. R., et al. 2005, MNRAS, 359, 1165 
Hayashi, M., et al. 2009, ApJ, 691, 140 

Hayashi, M., Shimasaku, K., Motohara, K., Yoshida, M., Okamura, 

S., & Kashikawa, N. 2007, ApJ, 660, 72 
Kashikawa, N., et al. 2004, PASJ, 56, 1011 
Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189 
Maiolino, R., et al. 2007, A&A, 472, L33 
Nakajima, T, et al. 2008, PASJ, 60, 435 
Okumura, S. K„ et al. 2000, PASJ, 52, 393 
Papovich, C, et al. 2007, ApJ, 668, 45 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., Scoville, N. Z., & Soifer, B. T. 1991, ApJ, 370, 158 
Schmidt, M. 1959, ApJ, 129, 243 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 
103 

Scoville, N. Z., & Good, J. C. 1989, ApJ, 339, 149 
Solomon, P. M., Downes, D., & Radford, S. J. E. 1992, ApJL, 398, 
L29 

Solomon, P. M., Downes, D., Radford, S. J. E., & Barrett, J. W. 1997, 
ApJ, 478, 144 

Solomon, P. M., & Vanden Bout, P. A. 2005, ARA&A, 43, 677 
Solomon, P. M., & Barrett, J. W. 1991, Dynamics of Galaxies and 

Their Molecular Cloud Distributions, 146, 235 
Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 

D. K., Reddy, N. A., & Hunt, M. P. 2004, ApJ, 604, 534 
Tacconi, L. J., et al. 2008, ApJ, 680, 246 

Tacconi, L. J., Genzel, R., Lutz, D., Rigopoulou, D., Baker, A. J., 

Iserlohe, C, & Tecza, M. 2002, ApJ, 580, 73 
Tsutsumi, T, Morita, K.-I., & Umeyama, S. 1997, Astronomical Data 

Analysis Software and Systems VI, 125, 50 
Young, J. S., Schloerb, F. P., Kenney, J. D., & Lord, S. D. 1986, ApJ, 

304, 443 

Young, J. S., Allen, L., Kenney, J. D. P., Lesser, A., & Rownd, B. 
1996, AJ, 112, 1903 



